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- - ■ Abstract 



The favoured R-parity violating SUSY scenarios for the anomalous HERA 
events correspond to top and charm squark production via the A' 131 and A' 121 
couplings. In both cases the corresponding electronic branching fractions of 
the squarks are expected to be <C 1. Consequently the canonical leptoquark 
signature is incapable of probing these scenarios at the Tevatron collider over 
most of the MSSM parameter space. We suggest alternative signatures for 
probing them at Tevatron, which seem to be viable over the entire range of 
MSSM parameters. 
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1. The Anomalous HERA Events 



The HI and ZEUS experiments have reported some anomalous high-Q 2 
events from the HERA e + p collider, which could be suggestive of new physics 
beyond the Standard Model (SM). The HI experiment has reportedly seen 12 
neutral-current events, e + p — > e + qX, at Q 2 > 15,000 GeV 2 against the SM 
prediction of 5 [1]; while ZEUS has reported 5 events at Q 2 > 20,000 GeV 2 
against the prediction of 2 [2]. Moreover the excess of 7 events observed by 
HI seem to cluster around a common e + q mass of 



which is not inconsistent with those of the ZEUS events [3]. Taken together, 
they represent an excess of 10 events, with a common mass range of 200-220 
GeV. They are based on the 1994-'96 data, corresponding to a combined 
luminosity of 34 pb' 1 for the two experiments, while the reported detection 
efficiency for each experiment is about 80%. This corresponds to a cross- 
section of 



for these anomalous events. 

These events have aroused a good deal of excitement in high energy 
physics in the past few months; and several possible mechanisms of new 
physics have been suggested [4-10]. It should be noted of course that the 
statistical significance of the signal is about 3a for each experiment, which is 
by no means conclusive. The ongoing experiments at HERA are expected to 
double the data sample in another year. While very welcome, this may not 
be sufficient to settle the issue conclusively. It is imperative therefore to see 
if the contending new physics mechanisms can be tested at other colliders 
- in particular at Tevatron, which has a considerable energy reach to probe 
these mechanisms. 

The new physics mechanisms suggested are mainly of two types — (i) a 
contact interaction term corresponding e.g. to a heavy Z' exchange [6,8,9], 
and (ii) the production and decay of a generic leptoquark, i.e. a hypothetical 
particle coupling to the e + q channel [4-8]. The first interpretation seems to 
be disfavoured on several grounds. The size of the contact interaction term 
required is at best marginally compatible with the upper limits from LEP and 
Tevatron colliders. Moreover it favours the standard e + q mass distribution, 



M ~ 200GeU, 



(1) 



a ~ OApb 



(2) 
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M = y/s ■ x, as given by the x distribution of the quark, instead of its 
clustering at a high value of M. For the same reason, the M integrated 
cross-section is expected to go down at large Q 2 . Instead, the clustering in 
M and the flat distribution over a very large range of Q 2 observed in the 
data clearly favour the formation and decay of a generic leptoquark in the 
e + q channel. 

Thus, it is natural to ask whether some of the extensions of the SM can 
naturally account for such a leptoquark. The leptons and quarks are unified 
in GUT, which naturally predict leptoquark states both as gauge (vector) 
bosons and Higgs scalars [10, 11]. However, the exchange of these objects 
in GUT generally leads to lepton and baryon number violating interactions, 
and in particular to proton decay. Thus the stability of proton implies these 
objects to be very heavy (M > 10 15 GeV), which puts them far beyond the 
reach of present or foreseeable future machines. While there are examples of 
GUT models like E$, having leptoquarks without baryon number violating 
couplings [10], there is no natural reason to expect them to be as light as a 
few hundred GeV. 

A more plausible candidate for a generic leptoquark in the mass range of a 
few hundred GeV is the scalar superpartner of quark (squark) in the R-parity 
violating SUSY model [12, 13]. In this case they can have lepton and baryon 
number violating Yukawa couplings and mediate proton decay as well. Unlike 
the gauge couplings, however, these Yukawa couplings are not constrained 
by any symmetry consideration. Therefore one can assume a finite value 
for the lepton number violating coupling while setting the baryon number 
violating ones to zero. The former ensures squark coupling to the e + q channel 
while the latter prevents proton decay. Thus in the R-parity violating SUSY 
model the squark can masquerade as a leptoquark and naturally account 
for the anomalous HERA events. Indeed these squarks are by far the most 
promising new physics candidates for the HERA events; and as such they 
have attracted a good deal of attention in the current literature on this 
subject [4-7]. The purpose of this work is to identify the most plausible R- 
parity violating SUSY scenarios for the anomalous HERA events and study 
the corresponding signals for the Tevatron collider. 

2. R-Parity Breaking SUSY Scenarios : 
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We shall consider the minimal supersymmetric extension of the standard 
model (MSSM) with explicit R-parity breaking [13]. The latter arises from 
the following Yukawa interaction terms in the Lagrangian: 



plus analogous terms from the permutation of the tilde, denoting the scalar 
superpartner. Here / and e (q and u, d) are the left-handed lepton doublet 
and antilepton singlet (quark doublet and antiquark singlet) fields; and i, j, 
k are the generation indices. The terms relevant for the HERA events are 



It is customary to assume a hierarchical structure for these Yukawa cou- 
plings in analogy with the standard Yukawa couplings of the quarks and 
leptons. The squark formation and decay processes corresponding to differ- 
ent choices of the leading A' coupling of eq.(4) are shown in Table - I. Only 
the 1st row corresponds to squark formation from a valence quark, while all 
other cases are from sea quarks. Knowing these quark fluxes one can easily 
calculate the cross-section for these processes at HERA for a given \'\/~B, 
where B denotes the squark branching fraction into the shown channel. The 
2nd column shows the sizes of corresponding X'VB, required to explain the 
cross-section (2) of the HERA events [5,6]. While the required size of the 
quantity is small for the valence quark case (1st row), it is larger by an order 
of magnitude in the other cases. Note that for the last two rows there is an 
equal probability of dj decay into the neutrino channel, i.e. B < 1/2. 

The last column shows the upper limits on these R violating couplings 
from other processes. The A' m limit comes from neutrinoless double beta 
decay [14], A' 112 113 limits from charge current universality [15], A' 121 131 limits 
from atomic parity violation [16], A' 122 133 limits from v e mass [17], A' 123 limit 
from forward-backward asymmetry [15], and A' 132 from Z-decay [18]. All but 
the last of these limits are taken from the recent compilation of ref. [5]. They 
are all la limits. 

More recently a precise measurement of atomic parity violation in 55 Ce 133 
has been reported[19]. The measured value of the weak charge is Q ex yy = 



— 72.11 ±0.27±0. 89, where the 2nd error is theoretical. This is in remarkable 
agreement with the SM value of Q SM W ~ -72.9 [20]. From 



L — XijkkljGk + KjJ-iQjdk + KjkdidjUk, 



(3) 



KjkihQjdk + kqjdk) + h.c. 



(4) 



Q w = -2[C lu (2Z + N) + C ld (Z + 2N)], 



(5) 
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1] 



it is clear that R-violating SUSY contribution can only add to the magnitude 
of Q SM w . Consequently a la bound would imply a very severe restriction 
on the coupling. But such a bound would not be justified since the data 
point is roughly la below the magnitude of SM. Therefore we have estimated 
the 2a bound, where the experimental and the theoretical errors of Q ex w have 
been added linearly. The resulting bound for a 200 GeV squark, A'lji < 0.10, 
is shown in the bracket in Table - I. This bound is consistent with the one 
obtain in ref.[21]. 

Comparing the 2nd and 3rd columns of this Table, we see that the re- 
quired coupling is reasonably small compared to its upper limit for the charm 
and top squark productions from a valence quark, i.e. 

e + d — > cl — > e + d, (7) 
e + d ^t L ^ e + d, (8) 

where the subscript L denotes left chirality. Moreover, the inclusion of NLO 
correction in eq.(7) and (8) is expected to reduce the required value of X'^y/B 
further by ~ 30% [22]. Therefore we shall concentrate on these two cases, in 
studying the R violating SUSY signal at the Tevatron. There is only one 
other case, where the required X'\/B is compatible with the upper limit on 
A'. This corresponds to top squark production from the strange quark sea, 

e + s ^t L ^ e + s. (9) 

This case has been recently studied in [18], where it was shown to give a 
consistent solution to the HERA anomaly over a limited range of the MSSM 
parameter space. We shall discuss this possibility while studying the top 
squark production scenarios of eq.(8) at the Tevatron collider. 

A brief discussion of the squark branching fraction B is in order. Under 
the assumption of hierarchical A' couplings, there is only one dominant R 
violating channel in a given scenario. The corresponding decay width of 
squark is 

r> =T I-A' 2 M. (10) 
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In addition the squark has R conserving decays into chargino and neutralino 
channels. The corresponding decay width Tr will be a function of the MSSM 
parameters, but independent of A'. Thus 

B = _Cl. (11) 

Note that the product \'\rB is constrained by the cross-section of the HERA 
events, as shown in the second column of Table-I. This can be used to elimi- 
nate A' from eqs.(10,ll), so that in any given scenario the branching fraction 
B is a unique function of the MSSM parameters. Thus 

where 

r^=-L(.04) 2 M (13) 

1 l07T 

for the favoured scenarios of eqs.(7,8). As we shall see below, in these cases 
one gets a B <C 1 over a large part of the MSSM parameter space. The signal 
for charm and top squark production at Tevatron will depend sensitively 
on this branching fraction. Note that for the scenario of eq.(9) one has to 
replace the factor .04 by 0.3 in eq.(13). Consequently the branching fraction 
B remains of the order 1 throughout the allowed MSSM parameter space in 
this scenario. 

Let us conclude this section with a brief discussion of the MSSM parame- 
ters, relevant for our analysis [23]. We shall assume a common gaugino and a 
common sfermion mass at the unification scale. Consequently the masses of 
the SU(3), SU(2) and U(l) gauginos (g, W and B) at the electroweak scale 
are related via their gauge couplings, i.e. 

M 3 = (g 2 /g 2 ) M 2 ~ 3.3M 2 , 

Mi = (5g ,2 /3g 2 ) M 2 ~ 0.5M 2 . (14) 

Thus there is only one independent gaugino mass M 2 . Of course the elec- 
troweak gauginos, W and B, will mix with the higgsinos, resulting in the 
physical chargino (Wip) and neutralino (-^1-4) states. Their masses and 
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compositions will depend on two more parameters - the higgsino mass pa- 
rameter (/i) and the ratio of the two higgs vacuum expectation values (tan (5). 
Finally, the right and left handed squarks of the first two generations are ex- 
pected to have roughly degenerate masses and so also the sleptons. These 
two masses are related via the renormalisation group equations which imply 
[12] 

M| ~ Mf + 0.85M 3 2 . (15) 
After QCD correction the physical (pole) mass of the gluino is [24] 

M~ g = (1 + 4.2a s /7r) M 3 = 1.15M 3 < M q . (16) 

Because of the large top quark mass, the top squarks are expected to have 
lower masses, with the hierarchy [25] 

M 1r < M h < Mg. (17) 

Moreover the large top quark mass also implies large mixing between ti and 
t R . This can further reduce the mass of the lighter physical state 

t ± =t L cos9 + t R sin9. (18) 

It also implies that t\ has significant left and right handed components. Thus 
it is a natural candidate for the anomalous HERA events. Indeed the possi- 
bility of this top squark production has been suggested for several years as a 
promising R-parity violating SUSY signal at HERA [25] . Therefore, we shall 
first investigate the implications of this scenario for the Tevatron collider. It 
should be noted that in this case the A' should be replaced by X'cosO in the 
second column of Table - I as well as in eq.(10). However the eqs. (12,13) 
remain unchanged. 

3. The Top Squark (Stop) Scenario at Tevatron : 

The dominant mechanism for stop production at Tevatron are the leading 
order QCD processes of quark-antiquark and gluon-gluon fusion [26] 

qq^hh , gg^tJi- (19) 

The R violating Yukawa coupling has negligible contribution to the produc- 
tion cross section and hence it is not considered here. Indeed the above 
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production processes hold equally well for leptoquark production. It was re- 
cently shown in [27] that these LO processes, combined with the LO structure 
functions, reproduce the NLO cross-section to within 15%. But combining 
them with the NLO structure functions underestimates the cross-section by 
a factor of 1.5 - 2. Therefore we shall use the LO structure functions CTEQ 
3L [28] in our analysis with the standard choice of the QCD scale Q equal to 
the squark mass. We have checked that one gets essentially identical results 
with the GRV 94 LO [29] structure functions. The structure functions are 
used via the PDFLIB version 6.06 [30]. 

Because of the large top mass, rrit — 175 GeV, the stop decay into the 
neutralino states can only proceed through higher order processes, 

ti -> cZi and t x -> bqq'Z 1 , (20) 

as long as 

M h <M Wi . (21) 

The corresponding decay widths are negligibly small [31] compared to the R 
violating decay width (13). Therefore, the dominant decay mode in this case 
is the R violating decay 

*i -> e + d. (22) 
On the other hand the R conserving decay 

h -> bWi (23) 

will dominate, when kinematically allowed. The corresponding decay width 
is given by [25] 



m^bW t ) = ^M h {l-MljMlf{Ai + A' R ) 



m b Ui 2 cos6 A Jr *___ n m t V* 2 smQ 



A R = ^-^ , Al = V*^ cos 9 H — 7=-^^ , (24) 

V2M W cos (3 V2M w s\nf3 V ' 

where U, V are the chargino mixing matrices and we have neglected the b 
mass in the phase space factor. 

Fig.l shows the resulting branching fraction for the R violating decay (22) 
over the relevant MSSM parameter space for = 200 GeV. The contours 
for the lighter chargino mass = 85 GeV and 180 GeV are also shown in 
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this figure. The former represents the limit of chargino mass which can be 
probed at LEP-II, although it has not been done so far for the R violating 
SUSY model. We see that the branching fraction for the direct leptonic 
decay (22) is B < 0.2 for a large part of the parameter space corresponding 
to M Wi < 180 GeV (going up to < 200 GeV for a 220 GeV stop). 

It should be noted that it corresponds to the most favoured region of the 
MSSM parameter space in terms of the naturalness criterion [32]. Thus the 
type of the stop signal at Tevatron will be sensitive to the choice of MSSM 
parameters. Let us analyse them case by case. 

(A) Direct Leptonic Decays : This corresponds to the direct leptonic 
decays (22) of the stop pair, resulting in a pair of hard and isolated e + e~ 
along with a pair of jets. Both the production and decay processes are 
identical to the leptoquark case recently investigated by the CDF and D0 
collaborations [33,34], except that in this case there can be no decay into 
the neutrino channel. Using a parton level Monte Carlo simulation we have 
estimated the signal cross-section following the CDF cuts [33]: 

(i) E Tel and E Te2 > 25 GeV, \i] e i\ or \i] e2 \ < 1, isolation (E^ c < 0.1E Te ) ; 

(ii) E T ji > 30 GeV, E Tj2 > 15 GeV, \r}ji, j2 \ < 2.5, cone size Rj = 0.7; 

(iii) M ee ^ 76 - 106 GeV; 

(iv) E Tel + E Te2 > 70 GeV, E Tjl + E Tj2 > 70 GeV. 

This is supplemented by the identification efficiency of 0.77 for the electron 
pair along with an azimuthal efficiency factor of 0.75 (corresponding to 85% 
azimuthal coverage for each electron). We have checked that the acceptance 
factor after each set of cuts agrees quite well with the CDF result. 

Fig. 2 shows the signal cross-sections before and after the above mentioned 
cuts against the stop mass. The right hand scale shows the number of events 
for the integrated luminosity of 110 pb~ l corresponding to the CDF data 
[33] . After these cuts there are 3 remaining events in this data against a SM 
background of similar size. The corresponding 95% CL limit of ~ 8 events 
implies a stop mass limit of ~ 200 GeV. One can eliminate the background 
by suitable cuts on the ej invariant masses, which can not be implemented 
however in our parton level MC simulations. The CDF group has achieved 
this while retaining 2/3rd of the signal by requiring the two ej invariant 



9 



masses to match within 20% and their average value to agree with the lep- 
toquark (stop) mass within a 3<r resolution error. This is illustrated by the 
dotted line in Fig. 2, which is obtained by multiplying the long dashed line by 
2/3 over the appropriate mass range. The corresponding 95% CL limits of 
3 events gives a better mass bound of M > 210 GeV. It may be added here 
that the D0 group has been able to eliminate the background while retaining 
80% of the signal via a hardness cut 

H T = E Tel + E Te2 + E Tjl + E Tj2 > 350GeV. (25) 

Combining this with an optimized set of kinematic cuts and a higher in- 
tegrated luminosity of 123 pb^ 1 they have reported a higher mass limit of 
M > 225 GeV [34]. 

Thus the CDF and D01eptoquark limits would rule out the stop scenario 
if the direct leptonic decays are dominant, i.e. > M^. It may be 

noted here that the alternative scenario of eq. (9) corresponds to B > 0.5, 
which may be already incompatible with the combined CDF and D0 data. 
However the favoured scenario (8) implies B < 0.2 over a large part of the 
MSSM parameter space corresponding to < 180 GeV (Fig.l). This 

means a reduction of the dilepton signal cross-section by a factor of at least 
25. Evidently the present data sample of CDF and DC 1 is in no position to 
probe a signal of this size. One expects a 20 fold increase in the integrated 
luminosity to 2fb~ 1 during the next (Main Injector (MI)) run of Tevatron. 
Moreover the increase of the CM energy from 1.8 to 2 TeV corresponds to 
a 50% increase in the cross-section. Thus one expects a 30 fold increase in 
the signal size during the MI run, which can probe the stop signal down 
to B — 0.2. Nonetheless there is a significant range of MSSM parameters, 
corresponding to B < 0.2, which will not be accessible to the dilepton signal 
during the MI run. With a further increase of luminosity to 20 fb~ l at TeV33 
it may be possible to probe the entire MSSM parameter space for the stop 
scenario via the dilepton channel. Nonetheless it is important to see if one 
can do better via the other decay modes of stop. 

(B) Mixed Mode : A more favourable signature for stop pair production 
in the low B( < 0.2) region is provided by the mixed mode, corresponding 
to direct leptonic decay (22) of one stop, while the other undergoes cascade 
decay 

fi -> bW u W 1 -> q'qZ 1} Z x ^ bud + h.c. (26) 
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Note that the R violating decay of the lightest superparticle (LSP) in this 
case is restricted to the neutrino channel due to the large top mass. Thus 
the final state will consist of a very hard e ± ) a large number of jets including 
a pair of b, and a modest amount of missing p T carried by the neutrino. We 
have estimated the signal cross-section for the following set of cuts: 

(i) p e T > 40 GeV, 1 77c | < 1, isolation (E$ c < O.lpf,); 

(ii) rij > 3 with E 3 T > 15 GeV, \rjj\ < 2.5, cone size Rj = 0.7; 

(iii) M T (e, E T ) ^ 50 — 110 GeV 

(iv) > 16-jet with E\ > 20 GeV and \r] b \ < 2. 

The 3rd cut is designed to reduce W decay background. We assume electron 
identification and 6-tagging efficiencies of 

e e . id = 0.8 e b = 0.3. (27) 

With the above acceptance cuts and 6-tagging, the dominant SM back- 
ground is expected to come from it production. Fig. 3 shows the signal cross- 
section for a 200 GeV stop along with the it background at y/s = 2 TeV. 
The MSSM parameters used are 

M 2 = 150GW, fi = -AOOGeV, tan/5 = 2 M w = 158GeV, M~ z = 77GeV; 

(28) 

but the signal cross-section is insensitive to these parameters. The signal 
shows a much harder electron px distribution than the it background. Besides 
one can exploit the clustering of invariant mass of the electron with the 
hardest jet at M ~ 200Gel / , to distinguish the signal from the background. 
Thus it seems feasible to separate the signal from the background while 
retaining the bulk of the signal size. 

The size of the signal cross-section in Fig. 3 is about 50 fb, corresponding 
to ~ 100 events at the MI run. But it is yet to be multiplied by the branch- 
ing factor ~ 2B. The smallest branching fraction over the allowed MSSM 
parameter space (Fig.l) is B ~ 7%. Hence this signature should be able to 
probe the stop signal at the MI run over the full parameter space of MSSM. 
An interesting feature of this signal is that the fe-jet pair will contain unlike 
as well as like sign fo's with equal probability, though it may be hard to test 
it experimentally. 
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(C) Cascade Decays : The largest event rate for B < 0.2 corresponds 
to the cascade decay (26) of each of the stop pair. The resulting final state 
consists of 4 6-quarks and a missing p T carried by the neutrinos; but the latter 
is severely degraded compared to the R conserving case. Consequently the 
process suffers from a large background from it as well as bb production. The 
leptonic decay of one of the charginos {W\ — > IvZi) will give a lepton (e,/i) 
40% of the time; but its detectibility will depend on the chargino-neutralino 
mass difference. In any case the it background remains and is 30-40 times 
larger at the level of the raw cross-section. Therefore one would need tripple 
b tag to separate this signal from the background. We have estimated the 
signal cross-section for the two channels with the following cuts: 

(I) j> T > 40 GeV, n b > 3 with E b T > 20 GeV, \r] b \ < 2; 

(II) v 1 t > 15 GeV, M < 1, E$ c < 0.lp l T , p\ T > 20 GeV, n b > 3 with 
E\ > 20 GeV, | ?76 1 < 2; 

assuming a 6-tagging efficiency e b = 0.3. Table-II shows the cross-sections 
for these two channels for a stop mass of 200 GeV and a CM energy of 2 TeV. 
The cross-sections are shown for several values of the M 2 and /x parameters 
with tan (3 = 2. The cross-sections are similar in size at tan/? = 10 as well. 
The first two rows correspond to the gaugino dominated region (M 2 <C |/x|), 
characterised by M% ~ M^jl. The last two rows correspond to the hig- 
gsino dominated region <C M 2 ), characterised by ~ M^. The first 
case implies harder lepton pt but relatively soft missing-j>r, compared to the 
second. Consequently the missing-pr signal goes up while the leptonic signal 
goes down as we go from the gaugino region to the higgsino one. Notwith- 
standing this complementarity, however, the absolute size of the signal is too 
small to provide a viable signature for an integrated luminosity of 2fb~ 1 . In 
short, we found no viable signature for detecting stop pair production, when 
both of them undergo cascade decay via (26). 

Note that the above conclusion is based on the current b tagging efficiency 
of 30%. An increase of this efficiency to 50%, as envisaged for the Main 
Injector run, will result in a 4-5 times increase in the signal cross-section. 
This could make it viable over a large range of MSSM parameters. 



4. The Charm Squark (Scharm!) Scenario at Tevatron 
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In this case one expects 8 species of roughly degenerate squarks along 
with a gluino of comparable or smaller mass. Only one of them, the left 
handed charm squark cl, has the R violating decay mode as required to 
explain the HERA anomaly 

c L -> e + d. (29) 

In order to estimate its branching fraction B, let us note that the width for the 
largest R conserving decays into the W dominated chargino and neutralino 

states are 

M\jMtfNl. (30) 

The masses and compositions of these states correspond to 

M M -M 2 <^M d and V n ,N i2 ~l. (31) 

Combining these with eqs.(12,13) one sees that the branching fraction for 
the R violating decay (29) is 

B < 1/20. (32) 

Therefore the direct leptonic decay channel (29) will not give a viable SUSY 
signature in this case. Instead one has to consider the cascade decay of the 
squarks and gluinos into the LSP. Fortunately the R violating decay of the 
LSP into e ± , 

~ A' 

Z 1 -H 1 ce + d(svd) + h.c. (33) 

provides a distinctive like sign dilepton (LSD) signature in this case. Using 
this signature one can test this R violating SUSY scenario over most of the 
MSSM parameter space even with the present Tevatron data. 

In this case one has to consider a host of production processes [26] 

w{gg) -»• m, w{gg) -»■ 99, qg(qg) -»• qg(^g). (34) 

For Mg < Mg, the cascade decay proceeds via the electroweak decays of 
squarks. Over most of the parameter space of interest the W± and Z 2 states 



Y(c L ^sWi) = —MAI 



r(c L - c zo = Jn M 5(i 
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are dominated by the W component while the Z\ is dominated by B. Thus 
the dominant decay modes of the left-handed squark are 

q L -> q'W u qZ 2 {W x -> qq'Z u Z 2 -> gg^), (35) 

while the right-handed squark decays directly into the LSP, 

q R -> gZi. (36) 

For < M<j, the cascade decay proceeds via the 3-body decays of gluino, 
i.e. 

g - ggWi(50%), gg^ 2 (30%), ^(20%), (37) 

where the approximate branching fractions are indicated in the bracket [35]. 
Finally the produced pair of LSP undergo the R violating decay (33). 

Thanks to the Majorana nature of the LSP (Zi), the final state di-electron 
have equal probability of having unlike and like signs. The latter constitutes 
a viable signature due to the low SM background in this channel. Indeed the 
LSD signature for R violating SUSY model has been considered in [35,36] for 
a variety of the A and A' couplings. We shall concentrate here on A' 121 as the 
leading R violating Yukawa coupling, as suggested by the charm squark sce- 
nario. We have estimated the LSD cross-section for the following kinematic 
cuts: 

p l T > 15GeV, \r]i\<l, Isolation (E$ c < 5GeV) . (38) 

The SM background for these cuts has been estimated to be only 2.4 fb at 
the CM energy of 1.8 TeV, coming mainly from WZ and ti channels[36]. 

Fig. 4 shows the LSD signal cross-section at a CM energy of 1.8 TeV for 
different choices of the MSSM parameters. In each case the signal is shown 
for a common squark mass M~ q = 210 GeV, with M g = 150-240 GeV. The 
upper limit of M g is suggested by eqs. (15,16) [37]. The contributions from 
the three finals states of eq.(34) to the signal are shown separately. They 
include the contributions from the small leptonic components in W±, Z 2 — > 
Z\ decays. But the bulk of the contribution corresponds to the dielectron 
coming from the LSP decays. We have not included any efficiency factor 
for electron identification. But it is clear from this figure that, even after 
making allowance for the identification efficiency, one expects to see at least 
half a dozen like sign dielectron events with the present Tevatron luminosity 
of llOpfe -1 . We conclude this section with the hope that this data will be 
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analysed soon to probe the R violating SUSY model, and in particular to 
test the charm squark scenario for the anomalous HERA events. Indeed 
it should be possible to do this with the dilepton data even without charge 
identification, where one can controll the SM background by choosing suitable 
kinematic cuts. 

5. Summary : 

The two favoured scenarios for the anomalous HERA events in the R- 
parity violating SUSY model are the production of top and charm squarks 
off the valence d quark via the Yukawa couplings A' 131 and A' 121 respectively. 
We have studied the prospects of testing these scenarios at the Tevatron col- 
lider assuming MSSM with common superparticle masses at the unification 
scale. In this case the size of the required A' coupling and the corresponding 
decay branching fraction B can be independently estimated as functions of 
the MSSM parameters. We find B <C 1 for the charm squark, while the 
same is also true for the stop over a large range of the MSSM parameters. 
Consequently the canonical leptoquark signature of dilepton plus dijets is of 
limited value in probing these scenarios at Tevatron. We suggest alternative 
signatures, which seem to be viable over the entire parameter space of in- 
terest. We have also considered the alternative scenario for stop production 
from a strange quark via the A' 132 coupling. In this case A' > 1/2; and this 
scenario may be already in conflict with the combined Tevatron data via the 
leptoquark signature. 

We gratefully acknowledge discussions with Drs. V. Barger, A.S. Belayev, 
D. Choudhury, S. Eno, S. Majumdar, N.K. Mondal, N. Parua, K. Sridhar 
and G. Wang. 
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Table— I: Different squark formation and decay processes are shown along 
with the size of the corresponding R violating Yukawa couplings, required 
to explain the HERA events. The last column shows the upper limits on 
these couplings from other processes. 



Process Reqd. X'y/B X' limit (M g ~ ~ 200 GeV) 



e+d 


— > Uj — > e+d 


0.04 


^'m 


< -004, A' 121jl31 < .13(.10) 


e+s 


— > Uj — > e+s 


0.3 


^112 


< .06, A' 122 < .09, A' 132 < 0.6 


e+b 


— > Uj — > e+b 


0.6 


^113 


< .06, A' 123 < .55, A' 133 < .003 



e+u -> 4 -> e+w 0.3 A' m < .004, A' 112113 < .06 

e+c -> 4 -> e+c 0.4 A' 121 < .13, A' 122 < .09, A' 123 < .55 



Table-II: Stop cross-section in the cascade decay channels with and 
without a lepton for = 200 GeV, = 2 TeV and the cuts described in 
the Text. The results are shown for several gaugino and higgsino mass 
parameters with tan/3 = 2. 



M 2 


n 




M~ Z1 






(GeV) 


(GeV) 


(GeV) 


(GeV) 


( I) i> T + > 36 (11)1+ t> T + > 36 


150 


-400 


158 


77 


3.3 


2 


100 


-800 


105 


51 


2 


1.4 


100 


-100 


101 


51 


3 


1.5 


200 


-100 


112 


89 


6 


1.2 


300 


-100 


111 


95 


8 


0.8 
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Figure Captions 

Fig.l : The stop branching fraction (eq. 11) is shown as a contour plot in the 
Mi-, ^ plane for the stop mass Mi = 200GeV, tanf3 = 2 and (a) 0$ = 0° 
(b) 0i = —45°. Contours for the lighter chargino mass of 180GeV(solid 
line) and 85GeV(dotted line) are shown. 

Fig. 2 : The stop pair-production cross section(No. of events for 110 pb~ l 
luminosity) shown against the stop mass at the Tevatron collider energy 
of 1.8 TeV. The solid line corresponds to the raw cross section, while 
the long and short dashed lines correspond to the dielectron+dijet cross 
sections following the CDF cuts as described in the text. 

Fig. 3 : The signal cross section(solid line) corresponding to the mixed mode(B) 
is shown for stop mass of 200 GeV along with the tt background(dashed 
line) against the px of the electron at 2 TeV. The MSSM parameters 
are M 2 = 150GeV, fi = -AOOGeV, and tan(3 = 2. 

Fig. 4 : The gluino-gluino (solid line), squark- ant isquark( short dashed) and 
squark-gluino(long dashed) contributions to LSD cross section shown 
as a function of gluino mass for (fj,,tan(3) = (a) -100 GeV, 2 (b) -300 
GeV, 2 (c)-300 GeV, 10. 
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